Astronomy & Astrophysics manuscript no. cyg 


February 5, 2008 


(DOI: will be inserted by hand later) 





Low frequency radio monitoring of CygnusX-1 

and CygnusX-3 

M. Pandey 1 , A. P. Rao 2 , G. G. Pooley 3 , P. Durouchoux 4 , R. K. Manchanda 5 , C. H. Ishwara-Chandra 2 



in 
o 
o 

(N 

a* 

00 

OS 



> 
o 
m 
in 
o> 
o 
in 
o 

^ 

Or 

6 
to 

in 
c3 



1 Dept. of Physics, Mumbai University, Mumbai - 400 098, India 

2 NCRA, TIFR, Post Bag 3, Ganeshkhind, Pune - 411 007, India 

3 Cavendish Laboratory, University of Cambridge, Madingley Road, Cambridge CB3 0HE, UK 

4 CNRS FRE 2591/CEA Saclay, DSM/DAPNIA/SAP, F-91191 Gif sur Yvette Cedex, France 

5 Tata Institute of Fundamental Research, Mumbai-400005, India 



Received ; accepted 

Abstract. We present results of monitoring observations of the micro-quasars CygnusX-1 and CygnusX-3 at 0.61 
and 1.28 GHz. The observations were performed with Giant Meter-wave Radio Telescope, GMRT, between 2003 
June and 2005 January. Variable, unresolved sources were found in both cases. CygX-1 was detected in about 
half of the observations, with a median flux density about 7 mjy at each frequency. The results show clearly that 
there is a break in the spectrum above 1.28 GHz. The variations in CygX-1 may be due to refractive interstellar 
scintillation. Cyg X-3 was detected in each observation, and varied by a factor of 4. For this source, models of the 
scintillation suggest a very long timescale (of the order of 40 yr at 1.28 GHz), and therefore we believe that the 
variations are intrinsic to the source. 



Key words, binaries : close - stars : individual : CygnusX-1, CygnusX-3 
continuum : stars 



ISM : jets and outflows - radio 



1. Introduction 

The X-ray sources CygnusX-1 and CygnusX-3 are accret- 
ing X-ray binary systems which have been classified as 
micro-quasars. The X-ray emission in such systems is be- 
lieved to originate in the inverse Comptonization of seed 
photons in the accretion disc; the sources exhibit several 
spectral states classified as low-hard, intermediate and 
high-soft state (Cui 1999). The spectral analysis of the ra- 
dio emission from these X-ray binaries with black hole can- 
didates (BHCs) suggests that the radio emission arises due 
to synchrotron emission from the high energy electrons 
emitted in the jets (Hjcllming & Johnston 1988; Hjcllming 
& Han 1995) and the radio emission is quenched when 
a source is in high-soft X-ray state (Corbel et al. 2000), 
and detected in both low-hard and intermediate X-ray 
state. High-resolution mapping of X-ray binaries, particu- 
larly the well-known object GRS 1915+105, has demon- 
strated the presence of relativistic outflows (hence the 
term 'micro-quasars'). Such motions have been mapped 
in Cyg X-3 (Mioduszewski et al 2001), and a jet-like fea- 
ture has been mapped in CygX-1 (Stirling et al 2001). 
There has been a variety of other studies, mostly at high 



radio frequencies (above 2 GHz). However, the low fre- 
quency characteristics of CygX-1 in particular have not 
been explored so far. 

To understand the physical mechanism connecting the in- 
flow and outflow of matter in such systems, a compre- 
hensive study of the correlation between radio and X- 
ray emission in different states of these sources is neces- 
sary. We have carried out systematic monitoring of micro- 
quasars using the Giant Meter-wave Radio Telescope, 
GMRT (Pandey et al. 2004) at low frequencies. In this 
paper, we present the observations of CygX-1 and CygX- 
3, which demonstrate variations in the flux densities at 
low radio frequencies for both these sources. The observed 
flux density variations can be ascribed either to intrinsic 
mechanisms specific to each source, i.e. a process occur- 
ring around the compact sources themselves, or to refrac- 
tive interstellar scintillation. The refractive scintillation 
can be described as the flux density variations caused by 
"focusing" and "defocusing" of electromagnetic waves by 
large scale plasma density irregularities in the interstellar 
medium (Rickett et al. 1984). 
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2. Observations and Analysis 

The radio snapshot observations (of duration 30 min) 
were carried out at 0.61 GHz and 1.28 GHz with a band- 
width of 16/32 MHz using the Giant Meter- wave Radio 
Telescope, GMRT. The flux density scale was set by ob- 
serving the primary calibrators 3C286 and 3C48. Phase 
calibrators were interleaved with 25 min scans on each of 
CygX-1 and CygX-3. The sampling time was 16 s. The 
data recorded from GMRT was converted into FITS files 
and analyzed using the Astronomical Image Processing 
System (AIPS). A self calibration on the data was used 
to correct for phase related errors and improve the image 
quality. Tables 1 and 2 summarize details of observations 
using GMRT at the two frequencies between June 2003 
and Jan 2005. CygX-3 was detected in each observation 
and CygX-1 during about one- half of them. Columns 1, 
2 and 3 in both the tables gives the dates of observations 
and the measured flux densities or upper limits. Column 4 
gives the corresponding rms noise in the field of the radio 
image. For the observations performed during telescope 
maintenance/test time the background noise was higher 
because of the reduced number of antennas available. For 
each field there are some 6 background sources detected. 
We record in Tables 1 and 2 the measured flux density of 
one from each field, at an angular distance of 15' (CygX- 
1) and 22' (CygX-3). The measured flux densities of these 
'control' sources are consistent with constant actual flux 
densities and the observed noise levels. This result gives 
us overall confidence in the reliability of the system. 
To supplement the data from GMRT, we have also used 
radio data on CygX-1 and CygX-3 taken with the Ryle 
Telescope at 15 GHz. The longitudes of the GMRT and 
Ryle Telescope differ by 74deg (5 hr), and so the obser- 
vations are usually not exactly simultaneous; in the worst 
cases the nearest 15-GHz observation is separated by 2 
days from the GMRT observation. The radio light curve 
at 15 GHz for CygX-1 was averaged for 10 min, that for 
Cyg X-3 for 5 min. The typical uncertainty is 2 mJy + 
3% in the flux scaling. In column 6 of the Tables we have 
listed the flux density at 15 GHz for the nearest observa- 
tion at that frequency. Column 7 gives the spectral indices 
between the GMRT observation and 15 GHz (in the sense 
S oc v a ). Column 8 gives the X-ray spectral state of the 
source as determined from the 2-12 keV light curves of the 
sources taken from the ASM archival data. 

In Fig. 1 and Fig. 2 we show the first ever radio images 
of CygX-1 and CygX-3 at 0.61 GHz. The AIPS subroutine 
jmf it shows that the images are consistent with point 
sources. 

3. Discussion 

3.1. CygX-1 

CygX-1 is bright X-ray source with a mass in excess of 
7Mq (Gies & Bolton 1986) and a distance of ^2.5 kpc. 
The companion is an 09.7 lab super-giant (V1357 Cyg, 
HDE 226868), showing a strong stellar wind (Persi et al. 
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Fig. 1. GMRT image of CygX-1. The contour levels are 
1.3mJy x (-1,1,1.4,2,4,8). 
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Fig. 2. GMRT image of Cyg X-3. The contour levels are 
2.3mJy x (-1,1,1.4,2,4,8) 



1980) and almost filling its Roche lobe (Bolton 1972). The 
observed radial velocity for the system v sin i is about 76 
km s _1 (Gies & Bolton 1982). The emission in X-rays 
and in the radio band vary over a wide range of time- 
scales, including modulation at the 5*?6 orbital period; 
see, for example, Brocksopp et al (1999) and Gleissner 
et al (2003). The X-ray emission is usually in the spectral 
state described as 'low-hard' (with a non-thermal com- 
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Table 2. Flux density of CygX-3 between 2003 Jun and 2005 Jan 
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ponent extending to a few hundred keV); the relatively 
rare 'high-soft' state is dominated by a strong low-energy 
component. The high-frequency radio emission tends to be 
suppressed during the high-soft states (Hjellming, Gibson 
& Owen 1975) (see Fig. 4). While the X-ray source is 
in the low hard state, the radio emission in the cm-wave 
band is persistent, but the magnitude does vary by a fac- 
tor of about 5, with very rare outbursts to higher flux 
densities. These variations are much less spectacular than 
those of some other members of this class such as Cyg X- 
3 and GRS 1915+105. The radio spectral index a, de- 
fined in the sense S cx v a , is close to zero in the range 2 
- 220 GHz (Fender et al 2000). Jet emission in CygX- 
1 has also been observed at milliarcsec resolution dur- 
ing the low-hard state at 8.4 GHz (Stirling et al. 2001). 
This jet is also not as dramatic (in terms of rapid vari- 



ability and relativistic outflows) as in the cases of jets 
observed in a number of other microquasars e.g. CygX- 
3 (Miouduszeski et al 2001), GRS 1915+105 (Dhawan 
et al. 2000), IE 1740.7-2942 (Mirabel et al. 1992), and 
GRS 1758-258 (Rodriguez et al. 1992). Observations by 
Gallo et al (2005), however, suggest that the jet may be 
inflating a large (5pc) shell-like structure around the X- 
ray source, implying that its contribution to the outflow 
of energy in this system is dominant. Fig. 3 shows the 
low-frequency flux density from GMRT plotted against 
the closest available 15- GHz flux density. Different mark- 
ers are used for the two GMRT frequencies. There is 
a clear correlation between the high- and low-frequency 
data, somewhat complicated by the upper limits on some 
of the GMRT flux densities, but nevertheless convincing. 
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Fig. 3. Flux -flux plot for CygX-1 showing the low- 
frequency GMRT flux densities (O 0.61 GHz, x 1.28 GHz) 
and the closest 15- GHz measurement. 

3.1.1. The radio and X-ray association of CygX-1 

The 15-GHz radio light curve and the RXTE/ASM X-ray 
count-rate and hardness ratio HR2 for CygX-1 from 2003 
Mar to 2005 Jan are shown in Fig. 4 (daily averages for 
both of the X-ray parameters are used). The X-ray state 
of the system changed several times over this interval. For 
the majority of the time, CygX-1 was in the low-hard 
state, but with three distinct periods of the high-soft state 
lasting 1-2 months each. The suppression of the 15-GHz 
flux density during the high-soft states is well-illustrated 
by these data. 

There was one exceptional flare recorded in the 15-GHz 
data, on 2004 Feb 20.27, where the flux density reached 
a maximum of 140 mJy in an event which lasted about 
30 min. Unfortunately there was no GMRT observation 
at this date. 

3.1.2. Spectral and temporal behaviour of CygX-1 

CygX-1 is a persistent but variable flat-spectrum source 
over the range 2 to 220 GHz (Fender et al. 2000). It is 
clear from table 1 and Fig. 3 that CygX-1 was detected 
at low frequency in only one of the 8 observations while 
it was in the high-soft X-ray state. Observations in the 
intermediate and low-hard X-ray states all resulted in de- 
tections, apart from two which both had relatively high 
noise levels. Therefore we conclude that the high-soft state 
suppresses the radio emission at low radio frequencies as 
it does at high frequencies. There is no clear distinction 
between the radio spectral index of the intermediate and 
low-hard states, being typically a = 0.3. While all of these 
spectral measurements are subject to the additional un- 
certainty from the lack of precisely simultaneous observa- 
tions, we also conclude that the radio spectrum rises from 



the GMRT observing frequencies to the regime where it 
is essentially flat (2 to 220 GHz). Part of the variation of 
radio emission from CygX-1 is associated with the X-ray 
state of the system, and is therefore intrinsic. A further 
part may be associated with propagation through the in- 
terstellar medium - in particular, refractive interstellar 
scintillation (RISS) may introduce amplitude variations 
(Rickett, 1990). We have used the NE2001 code made 
available by Cordes & Lazio (1993) to estimate the prop- 
agation parameters for CygX-1 and CygX-3. 

The NE2001 code estimates that the transition fre- 
quency (above which the scattering is weak, and there- 
fore the intensity fluctuations small) is around 6 GHz. The 
angular broadening is estimated as 0.3 mas at 0.61 GHz 
and 0.06 mas at 1.28 GHz. The scattering disc is therefore 
about 1 x 10 n m at 0.61 GHz and 2 x 10 10 m at 1.28 GHz. 
If we assume that typical transverse velocities involved 
(the source or the medium) are near 100 kms -1 , we derive 
timescales of 5 and 1 days at the two frequencies. These 
are clearly in the same area as those of the observed fluc- 
tuations, so it is quite possible that RISS is a contributing 
factor. We note, however, two further points. The angu- 
lar size of the (transient) radio jet mapped at 8.4 GHz by 
Stirling et al. (2001) is 10 mas: the low-frequency emission 
may originate in a region of this size (or larger); and the 
scintillation bandwidths estimated by the NE2001 code 
are 30 and 100 MHz at 0.61, 1.28 GHz. Both of these would 
reduce the observed RISS contribution to the fluctuations. 

3.1.3. Proposed geometry and emission mechanism 

From the above discussion, we conclude that the radio 
spectrum of CygX-1 consists of two parts, the persistent 
high frequency emitting component and the variable low 
frequency component. Both of these components imply 
the presence of a continuous jet with discrete plasmoids 
or clouds of electrons, emitting synchrotron photons. The 
observed behaviour of the source at radio frequencies can 
be explained by adiabatically expanding plasmoids in con- 
ical jets (Hjcllming et al. 1988), where the low frequency 
radio emitting component is significantly altered by re- 
fractive interstellar scintillation giving rise to high vari- 
ability in the radio flux density. A sketch of the proposed 
geometry is shown in Fig. 5. In our model compact plas- 
moids are ejected continuously at the base of the jet and 
expand freely outward. The low frequency emission from 
the newly ejected plasmoids will be absorbed due to syn- 
chrotron self absorption in an optically thick medium be- 
cause of its highly compact nature and high magnetic field. 
The plasmoid will be still visible at high radio frequencies 
(Miller- Jones et al. 2004). As the plasma expands under 
its internal pressure, it will become optically thin and emit 
over the whole radio-frequency band, although at a lower 
intensity. 

The observed flat radio spectrum of CygX-1 can thus 
be explained as due to superposition of individual self- 
absorbed synchrotron components emitted or plasmoids 
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Fig. 4. Cygnus X-1: Radio (15 GHz) (top) and RXTE/ASM (middle) light-curves; ASM hardness ratio HR2=(5- 
12)keV/(3-5)keV (bottom) during the interval MJD 52796 = 2003 Jun 06 to MJD 53392 = 2005 Jan 22. The dates 
of GMRT observations arc marked with vertical lines. 



with differing age profile and having longer life time. If a 
combination of size and the magnetic field are such that 
opacity effects are important only for the younger popu- 
lation, an inverted (or overall flat) spectrum is seen. The 
average spectrum can also mimic a flux level shift above 
a cut-off frequency, characterized by the optical thickness 
of the medium and the choice of age distribution. As dis- 
cussed earlier, the 15-GHz and GMRT flux densities are 
correlated; while we do not believe that the 15-GHz data 
are significantly affected by scintillation, this may not be 
the case at the GMRT frequencies and at least some of the 
variance in the 0.61 and 1.28 GHz data results form RISS. 
In conclusion, the radio emission in CygX-1 is mainly due 
to a continuous outflow of compact plasmoids or electron 
clouds into a jet and the spectrum consists of the compos- 
ite emission from each of these slow decaying ejecta and 
mini flaring events as observed in the source may result 
due to randomness of the age distribution of the ejected 
clouds. 
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Fig. 5. Model for CygX-1 jet 
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3.2. CygX-3 

Cyg X-3 is an X-ray binary system which does not seem to 
fit easily into the established classes of these systems. It is 
probably a low mass X-ray binary system. It has a 4.8-hr 
orbital period, and is located in the Galactic plane at a dis- 
tance of -10 kpc (Prcdchl ct al. 2002; Dickey 1983; Mason 
et al. 1979). The optical counterpart of the X-ray source is 
a WR star (van Kerkwijk et al. 1996), which is not visible 
in the optical band because of heavy interstellar extinction 
but is clearly seen above ~ 0.8^m. The optical data also 
shows a 4.8-hr modulation corresponding to the orbital 
period as seen in X-rays (Hanson et al. 1999). The source 
exhibits two spectral states low-hard and high-soft similar 
to CygX-1 and has been observed up to 500 keV during 
the flare mode. The quasi-periodic oscillations with peri- 
ods between 50 — 1500 s are the key characteristics of 
the source (van der Klis & Janson 1985). A detailed anal- 
ysis of the X-ray spectrum suggests that the total no of 
X-ray photons seems to be conserved at all times irrespec- 
tive of the state and the observed spectrum is consistent 
with a thermal source embedded in a hot plasma and en- 
veloped in a cold hydrogen shell (Manchanda 2002). The 
X-ray light curve of the source in the 2-12 keV band from 
the RXTE/ASM data shows frequent spectral changes be- 
tween the high-soft to low-hard states thereby suggesting 
large changes in the accretion rate on to the compact ob- 
ject. 

At radio wavelengths, Cyg X-3 is the most luminous 
X-ray binary in both its quiescent and flaring states 
(Waltman ct al. 1995). Huge radio outbursts have been 
reported in Cyg X-3 during which the flux density can in- 
crease up to levels of ~20 Jy; radio emission is suppressed 
( "quenched" ) to levels below 1 mJy for some days before 
large radio flares (Waltman et al. 1994). Jet-like structures 
with repeated relativistic ejection have been observed at 
various radio frequencies (e.g. Schalinski et al. 1998). On 
an arc-second scale, two-sided jets have been seen from the 
source in the N-S orientation, whereas a highly-relativistic 
(/3 > 0.81) one-sided jet with the same orientation has 
been reported on milli-arcsec scales with the VLBA (Marti 
et al. 2001; Mioduszewski et al. 2001). 

3.2.1. Temporal characteristics in the radio band and 
their association with the X-ray emission 

In the Table 2 above we have summarized the radio flux 
densities of the source as measured during various obser- 
vations with GMRT and from the Ryle Telescope data. 
Fig. 6 shows a flux-flux plot for Cyg X-3 for the GMRT 
and 15-GHz data, and Fig. 7 shows the 15-GHz and RXTE 
ASM data for the whole period. The timing of the GMRT 
observations is again marked with vertical lines in Fig. 
7. The last 4 GMRT observations, all at 0.61 GHz, were 
made during the high-soft state. The mean flux density at 
that frequency was 33.2 mJy, compared with 22.7 mJy for 
the 6 previous data points; the 15-GHz mean value was 
also higher, 150 mJy compared with 71mJy for the cor- 



responding 6 data points. We note that the radio/X-ray 
correlation is in the opposite sense to that for CygX-1. 
McCollough et al (1999) report both anti-correlations (in 
the quiescent state) and correlations (in the flaring state) 
between the hard X-ray flux (20-100 keV, as measured by 
BATSE) and the cm- wave radio flux density of CygX3. 

We investigated the possibility that RISS might be im- 
portant in the case of CygX3. The propagation conditions 
are more severe than for the case of CygX-1: the path 
length is much longer, and it has been known for some time 
(e.g. Wilkinson et al, 1994) that the scatter-broadening 
for this source is extreme. The NE2001 model is consis- 
tent: it suggests angular broadening of 20.59, 4.03arcsec 
at 0.61, 1.28 GHz. The corresponding timescales would be 
many years, and the narrow scintillation bandwidth (2 Hz) 
would suppress any observed scintillation. We conclude 
that RISS is not relevant to this study of Cyg X-3. 



Cygnus X-3 GMRT - 15 GHz 







X 










X 




X 




X 


o 


O 


o 


8 


o 






x 












O 


o 
o 











50 100 150 200 



S 15GHz/ mJ y 

Fig. 6. Flux-flux plot for Cyg X-3 showing the low- 
frequency GMRT flux densities (O 0.61 GHz, x 1.28 GHz) 
and the closest 15-GHz measurement 



To look for correlation between the radio emission from 
the source with its X-ray emission characteristics, we have 
plotted the RXTE/ASM X-ray light curve for Cyg X-3 in 
Fig. 7 along with the radio data. The timing of the GMRT 
observations is shown by the vertical lines. It can be seen 
that no large flares were observed during this interval. The 
radio emission is in the 'quiescent' state, typically 50 to 
200 mJy at 15 GHz. For the last month or so of this time- 
range, the X-ray spectrum softens: the RXTE ASM ratio 
HR2 falls consistently below 2, and the radio emission 
starts to become more erratic. This behaviour faded away 
after another month or so, and the source returned to the 
quiescent state. 
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Fig. 7. Cygnus X-3: Radio (15 GHz) (top) and RXTE/ASM (middle) light-curves; ASM hardness ratio (bottom), 
HR2=(5-12)keV/(3-5)keV during the interval MJD 52796 = 2003 Jun 06 to MJD 53392 = 2005 Jan 22. The dates of 
GMRT observations are marked with vertical lines. 



3.2.2. Spectral behaviour and the emission geometry 

As seen from the data in Table 2 and Fig. 6, CygX-3 is a 
persistent radio source at all wavelengths. CygX-3 is more 
luminous at higher frequencies. The data in Table 2 clearly 
indicates a low frequency turn-over in the source spec- 
trum below 1 GHz. As discussed earlier, such behaviour 
can arise due to synchrotron self absorption of the com- 
pact radio emitting plasma in an optically thick medium. 
The observed variability of the flux density is consistent 
with the assumption of a discrete ejection/plasmoid in 
adiabatic expansion. The ejection rate and/or lifetimes of 
plasmoids in CygX-3 are probably lower/shorter than in 
CygX-1 with little or no overlap of the spectra associated 
to the single discrete ejections. The spectrum and the fact 
that the radio emission is persistently at an high level may 
therefore imply an uniterrupted low-rate ejection. 

4. Summary 

In this paper we have presented the first ever low fre- 
quency maps of CygX-1 and CygX-3 along with the tem- 



poral behaviour from the monitoring observations of these 
sources over two years at various frequencies. It was noted 
that while Cyg X-3 is persistent source at low frequencies, 
CygX-1 is transient in nature with mean radio flux being 
7mJy and the source is detected at low radio frequen- 
cies only when the high frequency radio flux crosses the 
mean threshold of 15 mJy The spectral characteristics of 
both sources too have many similarities e.g. the spectral 
turn over at lower frequencies. Continuous blob of plas- 
moids arc emitted within the jet medium in both the cases. 
However the life time of the individual blob in the case of 
CygX-1 is larger than in the case of Cyg X-3 and hence a 
flatness in the radio spectra is seen in the case of CygX-1. 
The positive detection of CygX-1 at 0.61 to 1.2 GHz on 
several occasions along with the spectral turnover at low 
frequencies conclusively support the synchrotron origin of 
the radio photons contrary to its ambiguous thermal or 
non-thermal nature concluded from earlier measurements 
(Fender et al. 2000). In the case of Cyg X-3, a close corre- 
lation is seen in the change of X-ray hardness ratio and the 
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radio flux densities during the flaring and quiescent stages 
thus providing a consistent picture of disc-jet connection. 

Rigorous calculations were performed to look for the 
effect of refractive interstellar scintillation at low fre- 
quency on the data of these two sources. The present data 
clearly indicate that the source visibility above a threshold 
luminosity seen in CygX-1 is due to this effect. However 
the variability at high frequencies is intrinsic to the source 
and least affected by interstellar scintillations. Finally, we 
have discussed a common jet emission model for both the 
sources in which synchrotron radiation-emitting plasmoids 
originate at the base of the conical jet and the observed 
spectral features are due to the superposition of many such 
plasmoids with different age profile. Thus detailed spec- 
tral measurements in the m-cm regime (low frequency ra- 
dio) combined with simultaneous X-ray data can provide 
a clear understanding of the emission mechanisms in these 
sources. 
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